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Abstract 

A  linearized  Krook  equation  with  the  addition  of  a  point 
source  term  is  used  to  find  the  steady  flow  field  in  the 
vicinity  of  an  object  in  nearly-free  molecular  flight  through 
a  non-ionized  gas.   The  point  source  is  related  to  the  shape 
of  the  object  and  to  the  boundary  conditions  both  at  its  sur- 
face and  at  infinity. 

A  knowledge  of  the  flow  field  and  thence  of  the  distribu- 
tion function  near  the  body  permits  the  drag  to  be  determined 
as  a  first  order  perturbation  on  the  free-molecular  value. 

The  cases  of  a  plate  and  of  a  sphere  undergoing  specular 
and  diffuse  reflection,  respectively,  are  worked  out  in  detail. 
For  specular  reflection,  the  drag  on  a  plate  is  evaluated  for 
all  Mach  numbers;  for  diffuse,  due  to  the  length  of  the  calcu- 
lation, only  the  case  of  a  sphere  at  high  Mach  numbers  is 
considered.   The  results,  in  both  cases,  are  valid  to  first 
order  in  the  Knudsen  number « 
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NYO-10427 
Drag  on  a  Body  in  Nearly-Free  Molecular  Plow 

I.   Introduction 

The  region  of  nearly-free  molecular  gas  flow.  I.e.,  the 
transition  region  from  free  flow  to  the  start  of  Intermolecular 
collisions,  is  analyzed  in  this  report  with  a  view  to  computing 
the  drag  on  an  object  traversing  a  rarefied  atmosphere  at  a 
constant  speed.   The  method  to  be  presented  departs  from  those 
previously  used  for  this  problem,  namely,  the  "first  collision" 
and  the  "integral  equation"  methods  which  we  shall  now  describe 
briefly. 

In  the  "first  collision"  method  [1-7]  it  is  assumed  that 
the  most  important  correction  to  the  free  molecular  regime 
arises  from  collisions  between  free  stream  particles  that  have 
been  reflected  from  the  object  and  those  that  have  not.   Other 
types  of  collisions  are  disregarded  as  being  of  second  order. 
Baker  and  Charwat  [5]  have  analyzed  the  drag  on  a  sphere  using 
this  and  additional  assumptions,  which  are: 

1)  Billiard-ball  type  molecules. 

2)  The  molecules  reflected  from  the  surface  are  assigned  an 
average  speed  V  and  their  velocity  distribution  is  disregarded. 

5)  The  velocity  of  flight  V  is  sufficiently  large  for  the 
velocity  distribution  of  the  Incoming  molecules  to  be  disregarded, 

^)  V^o  «  ^  • 

5)  The  mean  free  path  of  the  emitted  beam  relative  to  the 


k   - 


Incoming  is  constant. 

6)  The  energy  accommodation  coefficient  is  unity  and  the 
body  is  cold. 

The  drag  coefficient  for  a  sphere  is 


V         V 
(I.l)     c^  -  2[1  +  0.44  ^  -    (0.46  =2.  +  1.3)apQSg^] 


where   a  is  the  radius  of  the  sphere,   p   the  number  density 

of  the  oncoming  stream,  and  S   the  collision  cross-section 
to         '  go 

of  the  emitted  stream  with  respect  to  the  oncoming  one. 

For  a  flat  plate  oriented  perpendicularly  to  the  flow. 
Baker  and  Charwat  obtain 


(1.2)     cj^  =  2[1  +  |!^-  0.32  ap^I^S^^  -  1.36  P^a  S^^] 

o  Ve 


Using  similar  assumptions  but  slightly  different  computa- 
tional methods,  Hammerllng  and  Kivel  [6]  arrive  at  the  following 
expression  for  the  energy  imparted  to  a  sphere  in  unit  time: 


(1.3) 


2  1       2 
^^   2  PoVo 


'-'\r 


(1  -0.2  =^  ap  S^^) 
y   ^o  eo^ 


Comparison  shows  that  0.44  V  /V   in  Baker  and  Charwat ' s 

c,  is  replaced  by  (l  -  V  /V  )   in  Hammerllng  and  Kivel 's,  and 

that   0,46  V  /V  ap  S    is  replaced  by  0.2  V  /V  ap  S   . 
o^e^oeo       ^       "^  o'^eoeo 

These  large  differences  demonstrate  the  sensitivity  of  the  results 
to  the  approximations  used  in  the  numerical  work.   For  example. 
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the  distribution  of  emitted  molecules  is  approximated  by- 
slight  ly  different  functions  in  each  case;  the  resulting  dif- 
ferences in  the  coefficients  is  due  to  the  fact  that  they 
represent  the  small  net  between  the  gain  and  loss  of  molecules 
hitting  the  sphere  due  to  collisions.   These  results,  although 
semi -quantitative ,  predict  trends  which  agree  with  experiments 
[7,8]. 

In  order  to  analyze  the  range  of  validity  of  some  of  the 
assumptions  made  by  Baker  and  Charwat ,  Willis  [9]  has  integrated 
the  Boltzmann  equation  formally  by  writing  it  as  a  pure  integral 
equation  and  then  solving  it  by  successive  approximations.   The 
series  thus  obtained  has  been  proved  convergent  for  the  special 
case  of  linearized  Couette  flow  [10];  furthermore,  it  has  always 
led  to  physically  reasonable  results.   By  approximating  the  . 
Boltzmann  collision  term  by  a  modified  KrOok  term  Willis  has 
shown  that : 

1)  The  collisions  between  streams  of  reflected  particles 
have  only  a  small  effect  on  the  drag  compared  with  the  first  col- 
lisions.  However,  this  effect  increases  as  V  /V   decreases. 

2)  Extrapolation  of  the  first  collision  results  to  higher 
values  of  V  /V   is  reasonable  for  pseudo-hard  sphere  molecules, 
but  is  not  justifiable  for  pseudo-Maxwellian  mqlecules.   Thus, 
the  choice  of  a  molecular  model  is  of  great  importance  if  the 
energy  accommodation  is  small  and,  therefore,  the  reflected  mole- 
cules are  no  longer  cold. 

3)  The  qualitative  agreement  between  the  two  methods  improves 
as  Vg/V^  ^0. 
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A  third  method,  as  yet  not  fully  analyzed,  has  been 
advanced  by  Enoch  [11].   It  is  applicable  to  high  speed  flows 
for  which  the  distribution  functions,  which  are  highly  peaked 
about  some  average  velocity,  may  be  approximated  by  6-functions. 
The  full  collision  integral  is  then  replaced  by  a  simpler  model 
which,  however,  retains  the  three  conservation  laws.   The  fore- 
going analysis  has  been  successfully  applied  to  high  speed 
plane  Couette  flow,  the  distribution  function  and  its  first 
three  moments  having  been  obtained  analytically. 

The  method  described  in  this  report  is  likewise  based  on 
replacing  the  Boltzmann  equation  by  a  simpler  model  equation 
which  yet  retains  its  essential  features.   The  equation  chosen 
is  the  "single-relaxation  time"  Boltzmann  equation  or  "Krook 
equation"  developed  by  Bhatnager,  Gross  and  Krook  [12]  and  which 
preserves  the  three  conservation  laws.   The  same  equation  was 
used  by  Grad  [13]  in  investigating  the  asymptotic  behavior  of 
the  flow  far  from  the  source.   In  the  asymptotic  case,  however, 
the  drag  could  not  be  found  since  the  source  function  was 
unknown.   A  knowledge  of  this  function  can  be  obtained  only  by 
determining  the  flow  in  the  vicinity  of  the  object  (less  than 
one  mean  free  path  away)  and  yet  at  a  sufficient  distance  for  the 
object  to  be  replaceable  by  a  point  source »   This  is  the  region 
to  be  studied  in  this  report. 

In  the  limit  of  a  point  source  and  granted  the  validity  of 
the  Krook  equation  we  obtain,  in  this  report,  the  first  order 
perturbation  on  the  free  flow.   This  result  is  rigorous,  i.e.,  it 
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encompasses  all  collisions  between  particles  and  not  merely 
"first  collisions." 

Our  results  are  obtained  by  following  the  method  used  by 
Grad  for  the  asymptotic  flow.   Briefly,  we  linearize  the  Krook 
equation  about  the  Maxwellian  distribution  at  infinity  and 
then  find  the  Fourier  transform  of  the  perturbed  distribution 
function.   This  yields  the  exact  flow  field  everywhere  in 
terms  of  the  Fourier  transforms.   The  inversion  is  different, 
however,  in  the  two  cases;  in  the  asymptotic  case  the  mean  free 
path  is  assumed  much  smaller  than  the  observer  distance,  where- 
as in  ours,  it  is  much  larger. 

These  methods  are  applied  to  finding  the  drag  on  a  sphere 
(and  concurrently  on  a  plate)  travelling  at  a  steady  speed 
through  a  gas  at  rest.   Calculations  for  two  different  boundary 
conditions  are  made,  namely,  specular  reflection  and  diffuse 
reflection.   In  both  cases  it  is  possible  in  principle  to  find 
the  drag  for  any  Mach  number,   However,  due  to  the  length  of 
the  computations,  only  the  following  cases  are  carried  through: 

a) .the  drag  on  a  plate  moving  in  a  direction  perpendicular 
to  its,  plane, for  all  Mach  numbers  with  specular  reflection  at  the 

boundary, 

b)  the  drag  on  a  sphere  at  high  Mach  numbers  only  with  dif- 
fuse reflection  at  the  boundary. 

We  obtain  the  following  corrections  to  the  drag: 
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D 


(1) 


4a 


u 


a 


'^    o 


— ^-^J^u^(.,  +  12k2)e   "+(l.,(l+uj 


u 


+  Skg] 


o 


2^ 


^    8/2 


dx  -  2u  (k,  +  8kp ) 7   for  the  plate, 

at  any  Mach  number 
undergoing  specular 
reflection. 


for  the  sphere  at  high  Mach  number, 
undergoing  diffuse  reflection, 


where 


u  =  Mach  number 
o 

A  =  cross-sectional  area  of  the  sphere  or  area  of  the  plate, 
a  =   radius  of  plate  or  sphere, 
L  =  mean  free  path. 


k^  -  -  16  /tt  , 

k  =  ^ 

7  =  ratio  of  temperature  of  Incoming  stream  to  reflected 
stream. 
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II.   Basic  Assumptions 

Our  aim  is  to  determine  the  drag  on  a  body  travelling  at 
constant  speed  u  through  a  rarefied  gas  at  rest. 

There  are  three  lengths  of  interest  in  this  problem:  the 
mean  free  path  L,  the  characteristic  body  dimension  a,   and 
the  observer  distance   R  from  the  object.   We  shall  assume 
that  R/a  »  1  and  that  L/R  »  1  which  together  imply  that 
L/a  »  1.   The  first  assumption  was  also  made  by  Grad  in  inves- 
tigating the  asymptotic  steady  flow  past  an  object;  thus  it  was 
possible  to  replace  the  body  by  a  point  source  function,  the 
value  of  the  drag  being  directly  related  to  this  function. 
However,  the  source  function  itself  could  not  be  determined 
because  Grad  assumed  L/R  «  1,  i.e.,  only  asymptotic  solutions 
were  sought  for  which  the  observer  was  many  mean  free  paths 
away  from  the  object.   However,  Grad ' s  results  are  valid  for 
any  value  of  the  Knudsen  number  L/a,  an  Interesting  conclusion. 
This  is  not  true  for  our  case. 

In  order  to  determine  the  source  function  and  thence  the 
drag,  we  must  know  the  distribution  function  f   in  the  neigh- 
borhood of  the  body,  i.e.,  within  a  mean  free  path  distance 
(l/R  »  l)o   However,  we  do  not  wish  to  go  too  close  to  the  body 
because  this  is  essentially  the  region  of  free-molecular  flow 
where  the  effect  of  collisions  is,  as  yet,  unimportant.   Now,  it 
is  precisely  the  effect  of  these  intra-molecular  collisions  on 
the  drag  that  we  wish  to  determine,  and  so  we  must  assume  that 
R/a  »  1. 
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This  last  assumption  allows  us  to  use  Grad ' s  procedure 
for  linearizing  the  Krook  equation,  and  we  obtain  the  same 
moment  equations.   However,  the  evaluation  of  the  moments  by 
Fourier  transform  techniques  and  the  subsequent  inversions  are 
radically  different  in  the  two  theories  because  of  the  different 
ranges  of  values  for  L/R.   The  linearization  procedure  and  the 
evaluation  of  the  Fourier  transforms  will  be  described  in  detail 
in  the  following  sections. 

In  order  to  treat  some  specific  problems  we  shall  consider 
the  cases  of  a  plate  moving  in  a  direction  perpendicular  to  its 
plane  and  of  a  sphere.   In  the  former  case  we  shall  assume 
specular  reflection  at  the  surface  and  in  the  latter,  diffuse. 
At  infinity,  the  distribution  function  will  be  assumed 
Maxwellian. 
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III»   Basic  Equations 

We  start  from  the  Krook  equation  with  a  point  source 
function,  namely. 

Here  f(x,^)  Is  the  distribution  function,   x  Is  the  position 

vector,   ^   the  velocity  vector,   v   the  collision  frequency, 

assumed  constant,  and  f   the  locally  Maxwelllan  distribution 

o 

function  given  by 

fo(^'^)  =  ^^7p  exp[-(^-u)^/2RT]  . 

f  (x,2,)  introduces  a  non-linear  term  into  equation  (lll.l) 
because  it  depends  on  f  through  p,  T,  u: 

p(x)  =  /  f  dl 
(III. 2)        pu  =  /  ^f  d^ 

p  =  p RT  =  /  i  (^-u)^f  di    . 

The  source  function  a (4)  will  be  subsequently  identified  with 
the  net  mass  flow  of  the  perturbed  distribution  function  through 
the  surface  of  the  body. 

Equation  (lll.l)  may  be  linearized  using  Grad's  method  [15] 
and  then  Fourier  transformed.   Although  repetitious,  this 
linearization  procedure  will  be  re-derived  here  since  it  has  a 
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fundamental  bearing  on  our  problem. 

We  shall  assume  that  the  distribution  function  f^"^  at 
infinity  is  Maxwellian;  the  corresponding  density,  velocity. 


.(o; 


and  temperature  are  denoted  by  p  ,  u  ,  T  .   Linearizing  about 
these  quantities  leads  to 


r 


(III-3; 


where 


f  ^  f^°^   +   g 


P    =   Po  + 

P 

u   =   u^  + 

u 

T   =   T^  + 

T    , 

.(o) 


(2TrRT^) 


^-372  exp[-(^-u^)72RT^] 


The  same  linearization  applied  to  f   yields 


(III.^) 


f 


(o){..^M 


c   _  3\  jr_  ,  c-u 

2RT     2''  T     RT 
o       00 


with   c  =  ^ -u 


Inserting  these  expressions  into  (lll.l)  results  in 


(111.5)   5  .  i  +  vg  =  vf  (°)[^  .  (jI^ 


)  ^  +  pm-  +  a5(x) 


2RT^    2^  T^    RT^ 
o       00 


which  is  linear  in  g   since 

P  =  /gd^ 


(III. 6) 


< 


PqU  +  P^o  =  /  ^  g  d4 

p  ^u  =  /  c  g  d^ 

P  =  /  3-  c  g  d^ 

p^RT  =  /  (3  c^-  RT^)g  d^  . 
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We   now   Introduce   the   following  dlmenslonless   quantities 
(except   for     L     and     x     where   L  =   l/v   Is   the  m.f.p.): 

P   =  p/pQ        '      T  =  T/T^        ,  p  =  p/p^ 

u  =  u//RT^   ,      u^  =  "c/^^ 

(III.  7)        I   =  i/f^   ,      c   =   c//RT^  =  !-u^ 

i  -    (RT,)^/VPo^      -  .  —1^  exp(-  I  c2)  =  (RT^)V2   ,0/^^ 

Then  equation    (ill. 5)    becomes 

(111. 8)  i.||+   1  i  .  g|p+  (1   c2-  |)T+c.u|+  a5(x)    , 
and  the   relations    (ill. 6)    take   the   form 

P   =  /i  dl    ,  u  =  /  eg  dl 

(111. 9) 

P  =  /  ^  c^'gdl    ,        T  =  /    (^  c^-   l)i  dl    . 

For  the  Fourier  transform  of  (ill. 8)  we  have 

(111. 10)  t(k,I)  =   ^^^]    P(k)+  (|  c^-  |)T(k)  +  c.u(k) 

1  +  1L|  •  k 


1+  lU  •  k 
where  a  function  /^(x)  and  its  transform  j2f(k)  are  Interrelated  by 
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1^{K 


e'^^'"^   /5(x)dx 


and 


>(x)  = 


(27r) 


3 


e^^"''  j25(k)dk 


-1  ~2 


If,  now,  we  multiply  (ill. 11)  by  1,0,  and  (^  c  -  l) 
and  Integrate  over  velocity  space,  we  obtain  a  set  of  five 
linear,  Inhomogeneous  equations  In  p ,  T,  and  u.   These  are 


(III. 11) 
/l  -  /fid! 


-/i(c^-3)nd!   -/ c  fid! 


\ 


-  /  i(c2-3)nd!  i-/^(c2-3)2nd|  -  /  ic^(c2-3)fid! 


-  /  c^fidl 


-/4c  (c  -3)ndl   5   -/c  c  fid 


2  r 


rs 


r  s 


T 


u 


=  L  a. 


(r) 


where 


/" 


n  = 


CD 


1  +  ILl-k 
a^(k)  = 


(III. 12)    < 


a2(k) 


a^^)(k) 


a(e)d! 
1  +  ILl-k 

((l/3)c^-  lja(l)dl 
1  +  IL^-k 


C-,  ad?, 


V, 


1  +  ILl-k 


The  subscripts   r  and   s   range  from  1   to  J>. 

The  next  step  Is  to  Integrate  the  various  components  of 
(ill. 11)  over  C,  solve  for  p,  T,  u  ,  and  then  perform  the 
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Fourier  Inversion.   If  L  Is  small,  the  denominator  of  Q   may 
be  expanded  In  powers  of  L  and  the  ^-Integration  and  subsequent 
Fourier  Inversion  of  (ill. 11)  presents  no  special  problems. 
This  Is  the  case  considered  by  Grad  and  subsequently  by 
Slrovlch  [1^];  the  observer  point  Is  many  mean  free  paths 
away  from  the  body.   In  the  near  region  In  which  we  are 
Interested,  however,  L  Is  large  and  the  | -Integration  of  the 
components  of  (ill. 11)  requires  special  consideration.   Some 
of  these  components  contain  a (I)  In  their  Integrands  and  It  Is, 
therefore,  appropriate  to  consider  the  nature  of  a (I).   This 
we  shall  do  In  Section  V. 
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IV".   Inversion  of  Equations  (ill.  11) 


The  5x5  matrix  (ill. 11)  Is  of  the  form 


(IV. 1) 


where 


l-£b.  .     eb.  . 
eb. .    1-eb .... 


1  =  1,2,. ..,5 
j  =   1,2,. ..,5 


and  e  is  a  small  parameter  (e^^  Vl)  . 

To  calculate  the  determinant  of  this  matrix  up  to  terms 

5 
in  e  would  be  extremely  laborious.   Fortunately,  due  to  the 

form  of  equations  (ill. 8)  it  will  only  be  necessary  to  calculate 

it  to  lowest  order  in  e.   For  completeness,  we  shall  compute  the 

determinant  of  (iV.l)  to  next  order  in  s.   When  solving  the 

linear  equations  (ill. 11)  for  p(k),  T(k),  and  u(k),  we  shall 

encounter  5x5  determinants  of  the  form 


x^ 


y^r 


X 


3 


X) 


£b^2 
1-eb 


^^13   ^^1^   ^^15 


22 


eb 


31 


1-eb 


33 


eb 


35 


More  precisely,  the  X;j_ '  s  may  occupy  any  column  but  the 
determinant  may  be  converted  to  the  form  shown  above  by 
changing  rows  and  columns. 
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In  order  to  solve  this  determinant  up  to  0(e),  It  Is 
sufficient  to  retain  the  terms  In  the  first  row,  the  first 
column,  and  the  diagonal.   All  other  terms  may  be  neglected. 
With  these  simplifications,  it  is  an  easy  matter  to  calculate 
p,  T,  and  u,  and  we  obtain 

(IV. 2) 

^  ^(^)  =a^^^k)(l+a^^)-a^-^a^(k)-a^g5g(k)-a^^a^(k)-a^^ap^(k) 


(2)  -  _        _  _ 

^  ^(^)  ^5^^^(k)(l+a^^)-a^^a^(k)-a^gag(k)-a^^a^^^(k)-a^^5p^(k) 

^  ^(^)^5p)(k)(l+a33)-a3^5^(k)-a3g5g(k)-a335^^)(k)-a3^5p)(k) 


We  note  that  all  the  a. .  and  all  the  source  terms  contain 
a  factor  1/l.   Therefore,  to  determine  the  moments  up  to  the 
necessary  order  as  functions  of  x,   it  will  be  sufficient 
merely  to  Invert  the  terms  containing  the  source  functions 
alone • 


-  18  - 


v.   Discussion  of  the  Source  Term 

For  L  — >  '=0,  the  source  term  o{i)    Is  precisely  the  free 
flow  perturbation  solution.   This  means  that  o{i)    at  any  point 
P(x)  Is  the  difference  between  the  stream  of  particles  reflected 
by  the  body  and  the  particles  that  are  unable  to  reach  P  due  to 

the  presence  of  the  body. 

2   2 
To  0(a  /R  ),  the  source  term  may  be  expressed  as 

(V.l)  o{l)    =   -      gl-dS 

where  the  Integral  Is  taken  over  the  surface  of  the  body,  the 
positive  normal  being  directed  Inward.   We  refer  to  the  article 
by  Slrovlch  [l4]  for  a  derivation  of  (V.l). 

It  Is  apparent  from  (V.l)  that  the  source  term  depends  on 
the  shape  of  the  body  and  on  the  boundary  conditions  on  g(x,I) 
both  at  the  body  and  at  infinity. 

As  a  first  approximation  we  assume  g,{x,i)    to  be  the  free 
flow  perturbation  solution. 

The  unperturbed  distribution  function  at  oo,  as  already 
mentioned,  is  Maxwell Ian: 


cd(|)  =  ^7p  exp  -  p  (I-u.) 


We  consider  the  following  cases 


■X- 

We  note  here  that  this  approximation  will  be  improved  upon 
by  using  an  iteration  procedure  which  we  shall  describe  in 
detail  in  a  later  section. 
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a)   Plate  -  specular  reflection 


|G< 


_  X: 


Figure  i 

For  a  plate  of  area  A  moving  in  a  direction  x-.    normal  to 
its  plane  we  have  (see  Fig.  1) 


(V.2) 


a(^)  = 
l?i|e 


A 


H(^l)+  Uile 


~      ~      ^      -W  Q 


H(-^l) 


where 


and 


H(^-j_)  =  1  when  ^-|_  >  0 


=  0  when  ?^  <  0 


In  expression  (V.2),  the  first  term  represents  the  distribu- 
tion of  particles  reflected  from  the  front  of  the  object; 
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from  this  we  must  subtract  the  distribution  of  particles 
Incident  on  the  back  which  Is  represented  by  the  second  term. 
The  first  two  terms,  therefore,  account  for  the  perturbed 
distribution  at  the  front.   Similarly,  the  next  two  terms 
represent  the  perturbed  distribution  at  the  back. 

b)   Plate  -  diffuse  reflection 
The  source  term  is 

(V.3)     a(i)  = 

where   7  =  T./T  ,  the  ratio  of  the  temperatures  of  the  incident 
and  reflected  streams,  and 

(V.4)  N^°)= 

2        [    /—I  -|u   cos  0       /— 1  _ 

^^X7^^  y?  u   cos  0  +  e        °  +   y^  u   cose  erf  (iu   cos   0) 

(2Tr)^/         /2      o  /2      o  ^-^    o  ^ 

b 

2  (         r-<  -iu^COS^0         /— ^ 

^172  1     /?  "0^°^  ^  ^   ^        °  +  yj  u^cos  0  erf  (iu^cos   0)| 


with 


erf (iu    )    =  -^ 

O^  v^ 


*^o 
/-      o        2 


e        dt 
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c )   Sphere  -  specular  reflection 


Figure  2 

It  is  easily  shown  that  the  source  term  is  given  by 


where  a  =  radius  of  the  sphere, 


d)   Sphere  -  diffuse  reflection 
(V.6) 


oil)    =  Tra^l  N^°)|ije-2'>'^\(e-^) 


(2Tr) 
1 


372  I'-i 


eje 


-i[(ei+u^)^+i^] 


i„,c2 


(°)|?  i.-W 


H(ei)+N^"^|^l|e-^1"=  H(-^^) 


(27r) 


572  l^il^  H(-^i) 
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It  Is  easy  to  check  that  In  all  cases  the  condition 

J  o{l)dl   =  0 
is  valid,  i.e.,  conservation  of  number  is  satisfied. 


VI.   Evaluation  of  the  Coefficients  of  the  Matrix  (ill. 11). 

A  knowledge  of  the  source  term  o{i),    as  discussed  in  the 

preceding  section,  enables  us  to  calculate  the  terms   a-,(k), 

ap(k),  cj^(k)  of  expression  (ill.  11).   After  evaluation  of  the 

a. .,  the  matrix  is  inverted  and  the  expressions  obtained  for 
_ij   _        _  ^ 

p(k),  T(k)  and  u(k),  Fourier-transformed  (Section  VII ).   In 
Section  VIII  we  shall  solve  equation  (ill. 8)  (i.e.,  the  linearized 
Krook  equation  for  the  perturbed  distribution  function)  where  the 
inhomogeneous  term  is  now  a  known  function  of  x. 

Although  we  shall  not  need  the  a. .(k)  as  such  for  our  compu- 
tation of  the  drag,  we  shall  yet  list  them  in  Appendix  I  where 
they  will  serve  for  calculating  the  source  terms. 

We  shall  now  describe  how  to  evaluate  a  typical  a.  .(k)  by 
integrating  over  C-   The  source  terms  may  be  integrated  over  %    in 
similar  fashion  as  they  are  merely  combinations  of  the  components 
of  the  a.  ..   We  have,  for  example 


(VI. 1)   a^^ 


fid!  =  ^7^ 


'^exp-i[(|  -u  )2  +  i2^L]   ~ 

^—^ ^  d^  , 

1  +  ILl-k 
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assuming  that  the  free  stream  velocity  at  oo  is  In  the  x-. 


direction. 

First,  the  k  coordinate  system  Is  rotated  so  that  \t 
only  has  an  x,  component;  hence  (c,  +u  )-k=c, Ikl  +u  -k 

1"^  v]_0'         1''      O' 

the  primes  referring  .to  the  new  coordinate  system  (in  the 
following  calculation  we  shall  drop  the  primes).   Then 


(VI. 2)  a^^   =  -^ 


1 


y2ir  lL|k| 


1-2 


dc- 


1  "o'^  -^ 

(c.  +  .•  T  I  V  I  +  1171—) 


'1      lL|k|       |k 
Expression    (V.2)    Is   of  the   form 

2 


(VI. 3)  Iq(z)    = 


e        dt 
(t-z) 


where      z   =  x  +   ly. 


which,  as  we  shall  show,  may  be  reduced  to  a  function  of  a 
tabulated  Integral  [15]- 


-z 


Let  t  -  z  =  -u,  and  assume  y  >  0;  then  I„(z)  -   e   f (z) 


where 


f(z) 


^    2   „ 
-u   -2zu 
e   e 


u 


du  , 


L 


the  path  of  Integration  L  being  shown  In  Figure  3. 


r 


Z=x+iy 


w 


L 


-00 


Figures 
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■^ 


f 


+  00 


-*X 


Now 

(VI. 

^) 

and 

(VI. 

5) 

f  (z)  =  -  2/¥  e' 


f(z)  =  -2/¥ 


e   dt  +  C  , 


o 


where   C   Is  a  constant  to  be  determined.   But 


00  „ 

r     -2 


llm       Io(y)   =     lirfi 
y  — >  °°  y — >  °° 


e        dt 
t-ly 


0    . 


Therefore,  since  Iq(z)  — >  0  as  y  — >  oo,   f(2)  =  lQ(y)e~^  when 
y  — >  00  must  also  approach  zero. 
It  follows  that 


(VI.  6)    C  =   llm  2/¥ 
y  ^~ 


yi 

e   dt  =  2)/¥  1 


-u  -, 
e   du  =  TTi  , 


(for  y  <  0  we  would  obtain  -ttI).   Therefore, 


(VI. 7)     Io^^(2)  =  e"^  IttI  -  2)/¥ 


z 

^  t2 

e""  dt  >   y  >  0 


(VI.  8)     I^^)(z)  =  e  ^  (-TTl  -  2/¥  1  e^^dt  ^  y  <  0 


o 


^  t2 


The  exponential  Integral  /  e   dt   is  tabulated  in  reference  15 . 

o 
In  terms  of  our  notation,  we  have 
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(VI.  9) 


^11 


/¥  lL|k| 


e  "  <|  TTl  -  2  /¥  1(a) 


where 


a  = 


/2  \L|k| 


Vl 
Ikl 


a   ^2 


with  a  positive  Imaginary  part  l/L|k|   and  l(a)  =  /  e  dt. 

o 
The  remaining  coefficients  a.  .  of  the  matrix,  together  with 

an  outline  of  their  evaluation,  are  given  in  Appendix  I. 

For  specular  reflection,  the  source  terms  (III. 11) 

Integrate  to 


k,A  2\  2 

(VI. 10)        a^(k)    =  — ^^—^  e  "   ^1/?  a  +  e"    -  2al(a) 


IL  k 


"   A  „2f  )        k^A  o2(  o2  , 

-^ e-°^  <^7ri-2/¥  1(a) ^  e'^  \l/¥  ^  +  e^  -2pi(p) 

/i  iL|k|  [  J      lL|kr 

^-^  e-n7rl-2/?  iO)      ^  a(^)(a,k)    -  a|i)(p,k) 


~fll               Ae""^      r(k2  +  k|))/7r  1 
(VI. 11)        5^1^(k)    ^^^^-^^       2        5 


IL  k 


+  /it  1  k^a^  +  k^ae°' 


2  .  ,  2 


-  21(a) 


(kp  +  k^)  2   p' 

"^g     ^     +   k^a^ 


Au  k,  2\  2 

+        °   ^   e-°^     l/¥  a  +    e°^    -  2Ia(a); 


iL|k| 
-  21  0) 


A  ^_p2|    (k^  +  k|)  >/¥l 


iL|k| 
+  /^  1  k^p^  +  k^pe^ 


(k^  +  k?)  p   pil      Au  k.         _^2j  g^ 

—2 2L  +  k?32    J  +  _^J_  e   P    l^l/rp  +  e^    -  2pi(p) 

ILlkr 


.^l)(^)(a,k)    -a^l)(i)(p,k) 
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(VI. 12)         o[^\k)    =  — ^2      g-a    |_  )^^  +  yr^  i^2_j_^ga  _^    ^^^^^^_2^2^ 
^  iL|kr 

+      ^     ^   p  1/Fa  +  e"    -  2al(a)      -  — L2      ^-p      _  ]^ 

ILlkl  (  J         ILlkp  ' 

P  o2  P    ")        u   Akpe"^  (  „2 

+  /¥  ip^+PeP  +   I(p)(l-2p2)      +      °     2  )  ^  y/^  p  +   eP 

ILlkl 


-2PI0)      .a(2)(r^)(,,,k)    -  a(2)(l)(p,k) 


a.^    '^  (k)    Is   the    same   as   o^      (k)   with  kp   replaced  by  k^. 


-  k  e"*^  A  (  2 

(VI.  13)         ap(k)    =  ^  — ^  ^  -  2/^  la  +  /?  la^  +  e°^    (a^-  |) 

"^  ^      iL|kp   (  "^ 

P    /         u        -a    .   (  P  2 

+  2al  (a )  (2  -a    )      +  — -  2  j/?  i  +  )/^  la     +  ae^ 

3    lL|k| 

+   2(2-a2)l(a)  }   -  ^  k     ^^-^  [  -  2 /If  ip   +  /?  ip^  +  e^^P^,  |) 

+  2piO)(2-p2)(   +  ^  el^A:l-  2)/T  1  +/¥  ip2   +  pe^ 

3    lL|k| 

+  2(2-p2)i(p)|   ^  S^(a,k)    -  a^(p,k) 


where 


-,     ,       .  u   k-, 

a   -     ^    ^      ^  °   1 


/2"  I  L|k|  |k| 

-,     f      .  u   k, 

/2      L|k|  |k| 
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A  may  represent  the  area  of  the  plate  or  the  equatorial 
cross-section  of  the  sphere.   The  meaning  of  the  functions 
a    (a,k)  and  a    (P,k)  Is  obvious;  the  superscripts   1   and  r 
refer  to  the  Impinging  and  reflected  streams  respectively. 
For  diffuse  reflection  these  source  functions  become 


(VI. 1^)        a^(k) 


(o)         (o)      ^1^      eT/2L2|k|2 

1^  iJ_j 


-7/2L^|k|^ 


Yyir     _  K^Y  1   J 
2L|k|         L|k|  \  )/2    L|k| 


-  5J^)o,k) 


(VI. 15)        S^^Mk)    =  A(n).°Un^°^)   ^- —   ^       2        ^ 


Yv  1  yk-, 
2iFhd^ 


^  I2  TlkT     1 


(k^+4) 


yk- 


L     k     • 


Y 


2    L  k 


~(1)(1) 
5 


0,k) 


(VI. 16) 


3<2)(lc)    = 


ILlkp         ^  "  [2         2L2|kr 


+ 


.1^^-/^^'^'^'- 


2  LTkl 


-a(2'<^)(P,k) 


(VI. 17) 


by  k^ 


Is  the  same  as  a^  '^  (k)  with  k^  replaced 
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,7/2L|k|^ 


'^  ^J   1     D        IL  kr    Lk   L''  k r 


-y/2L  k 


;i^^|-fr^^/?r^)(-il^j 


a^^^O^k)  . 


For  the  purpose  of  calculating  the  drag,  we  shall  only 
need  to  know  these  source  terms  up  to  0(1/L).   For  specular 
reflection  there  Is  a  considerable  simplification  in  the  source 
terms  if  this  approximation  is  used  and  this  case  will  be 
carried  through  for  the  plate.   For  diffuse  reflection,  however, 
the  source  terms  are  considerably  more  complicated. 

Simplifications  occur,  however,  if  we  let  u  — ^   co,  and 
results  for  a  sphere  undergoing  diffuse  reflection  will  be 
presented  in  this  report. 
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VII .   Evaluation  of  p(x)  for  Specular  Reflection 

In  this  section,  the  Fourier  Inversion  of  p(k)  will  be 
performed.   We  shall  only  consider  specular  reflection, 
reserving  diffuse  reflection  for  a  subsequent  section.   The 
calculations  will  be  performed  using  the  free  flow  source 
terms  discussed  in  Section  V.   The  results  obtained  here 
will  be  valid  for  any  Mach  number  and  up  to  0(i/l)  in  the 
Knudsen  number. 

We  shall  discuss  the  Fourier  inversion  of  p(k)  in  some 
detail.   The  calculations  for  the  other  moments  are  essentially 
the  same . 

A  typical  term  to  be  evaluated  is 


r 


(^"•1)        Poo^-)    -J^ 


e^^'""  5'^(k)dk  = 


SttXtt  Au^ 
(27r)^L 


u. 


k. 


2°    k?  4-1^2 


^W'F)^ 


1  - 


/2     (k2  + 


COS   k-^x-j^J^(;^R)lf 


where      R 


-i-- 


x^    +   X- 


,  ^=)^ 


2 
function  of  zeroth  order.   The  Integrand  of  (VII. 1)  is  now 


2  +  k^  ,   and  Jq(x)  is  the  Bessel 


expanded  in  a  series  which  may  be  shown  to  fulfill  Abel's  test 
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for  uniform  convergence  [l6].    Thus  we  may  Integrate  the 
series  term  by  term  over  k-,  . 

Intergratlon  over  k.,  yields 


(VII. 2)    p^Jz) 


Stt^Au. 


{2ir)\ 


(-1) 


2n+l 


n=o  nI2 


^2n(^^l) 
r(2n+i) 


^ 


2)/2     r(2n  +  |)   2n+l^    1^ 


using  expression  (11)  of  Appendix  IV  and  properties  (4),  (7)^ 
and  (8)  of  the  modified  Bessel  functions  K  (k),  listed  In 
Appendix  III . 

In  order  to  Integrate  over  k    we  must  again  show  that  the 
Integrand  Is  a  uniformly  convergent  series  In  ^ ,   Referring 
once  more  to  Abel's  test  we  have: 


* 


A  discussion  of  Abel's  (or  Hardy's)  test  for  uniform  conver- 


gence is  contained  in  [l6].   It  states  essentially  that  if  v  (x) 
be  either  monotonlc  increasing  in  n  for  each  fixed  x  in  (a,b) 
or  monotonlc  decreasing  in  n  for  each  fixed  x  in  (a,b),  then 
a  (x)v  (x)  is  uniformly  convergent  in  (a,b)  provided  that 


n 


n 


n 


1) 


a  (x)  is  uniformly  convergent  in  (a,b). 


n 


n 


2)   3  1^  such  that  |v  (x)  |  <  k  for  all  n  when  a  <  x  <  b, 
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~2  vH 


.)   I^^ 


(-1) 


3n+l 


n  -     ,-,     2^"+^  nan  +  |)   '"' 

are  uniformly  convergent  series  Independent  of   x,  . 

b)  The  K  (^x,  )  are  monotohlcally  Increasing  In  n  for  a 
given  R  X-,  . 

c)  Tf^""^"*"  ^2^(7^  x^)  and  4.^^'^^   ^2n+l^^^l^  approach  zero  as  4 
approaches  zero  from  properties  (l)  and  (2)  of  Appendix  III. 

d)  11m  K  (;qx-|)'^  e~   which  approaches  zero  faster  than  'n 
raised  to  any  power. 

We  conclude  that  all  the  requirements  of  the  test  are 
satisfied  and,  therefore,  that  we  may  Integrate  (VII. 2)  term 
by  term.   Using  expression  (12)  of  Appendix  IV  we  obtain  (up 
to  0(1/L)) 


(VII. 3)   Poo^^^  = 


^-^  Y     i^^[^fU^^F(n.l,l;l;4) 


K2         r(n+l) 


2 

X-, 


where  R  =  )/x|  +  x|  and  x^  =  r  cos  0  (see  Figure  2)  with 
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The  function  P(ajb;  c;  z)  is  the  hypergeometric  function, 
some  of  whose  properties  are  listed  in  Appendix  V. 

The  evaluation  of  the  other  terms  in  p(x)  is  entirely 
similar.   For  example,  we  have 


(VII. 4)        p..(x)    =  -^ 


r 


(27r)^     ^ —  nl  r(n+l) 


Xn 


^1 

n  r{n+2)  2      ^2        ^U2/ 

(2P+l)pI       r(n+p  +  |)  x^ 

y^f=^\y^i       (2EH-i)p:    r(M-p+|)  x^ 

where  A  may  refer  either  to  the  area  of  the  plate  or  to  the 
cross-sectional  area  of  the  sphere. 
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VIII.   Evaluation  of  p(x),  T(x),  and  u(x)  at  High  Mach  Number 
and  for  Diffuse  Reflection 

The  values  at  the  a.  .  and  the  source  terms,  for  high  Mach 
number,  are  listed  in  Appendix  II.  For  our  purposes,  we  shall 
only  need  the  source  terms  up  to  0(1/L).   These  are 

iL  k 


(VIII. 1) 


,U) 


AN 


(o)  ^_ 
f  /ir 


Au. 


1    ^l 


04-  (k)  =  -^ '-^  (k^  +  k^)  +  -—f^  k,  (1  - -|  ^  ) 


L  k 


^  Lkr 


iL  k 


^^0^2 
iL  k 


X2  |kl 
2 


1   ^1 


Y2     |k 


2 


ap)(k)  = 


A% 


'"   2 

—    /  X  k-,  k^    Au  k^        -,   k-, 

2    ^    L|kp    iLlk|2      /2  |k|2 

3  \ 


/,  \  _  c  I   1   1_    i_        1 

""^    ~  iL|k|  b)^  |kl  "  6  |kp 


Applying  the  Fourier  Inversion  and  following  the  methods  of 
Section  VII,  we  obtain 


(VI1I.2)   p^^(x).-^— ^ 


Ik.x  ~ 


CN 


(o) 


l^iTp/l  n  . 


R 


^(k)dk  =  ^-  r(f)F(f.i;i5 -%) 


J 


^1        ^  ^1 


r 


e^^""  o^{k)dk  = 


(VIII. 3)   T^^(x)  =  -^-^ 


3x' 


1 


j/2 


'2^  ^2 


xl       2  r(2) 
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(VIII. 4; 


<Vi^) 


i^ir) 


3 


e     oX    '  (k)dk 


.  C   ~  p(|),(3  ,^,^   R  )  ^  ^  £^[2^  ^^^  R 


x^     X2  r(2)     2 

^1  ^ 


Xn 


(VIII. 5) 


^<^o>(^) 


r 


L 


(27r) 


3 


Ik'X  ~(2  )  /,  \ -I, 
e     a^  ' (k)dk 


an|°V 


i::i2i 


P  — cos  ^  \   P(2,2;2;  -\  ) 

(27r)2  r(y2)       x^  x^ 


2Au 


(27r) 


2  /2 


°   Y^  cos 


r(2)r(2)-    xj      2 


R 


) 


+  2(1  -  -^)Au  ^^  cos/^  rd)  Af(2J;2;-^) 
y2    °  (27r) 


,-'"'0    ir.    n2  --"  a^  I  K^i       jxv^,2'^'    2 

^1  ^1 


(VIII. 6) 


OO   ^    '  /r^_\l> 


(27r)- 


e"^  "^  a"4   (k)dk   which  Is  the 


-(2) 


same  as  u'^  '^  (x)   with   cos  j6     replaced  by   sin  j6 , 


where 


4-7r  Xtt  a 


(27r) 


3 


TVT  (  O  )         7     ~ 

NX  '    -   7^—  u  cos 
I      dir     o 


0   when  u 
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IX.   Evaluation  of  i(x),  the  Perturbed  Distribution  Function  for 
a  High  Mach  Number  Sphere  Undergoing  Diffuse  Reflection 

We  start  from  equation  (ill. 8),  namely, 

!-||+  ^g  =  r(p(^)  +  (l^^-l)  •^^'^)  +  c.u(x)   +  5(!)5(x) 

of  which  the  right-hand  side  is  now  a  known  function  of  x. 
We  shall  need  to  know  p(x),  T(x),  u(x)  only  to  lowest  order, 
i.e.,  we  shall  use  p  ^(x),  T   (x),  u^^(x)  in  which  the  terms 
in  1/L  are  neglected. 

For  o(l)   we  use  expression  (V.6).   Transforming  to 
spherical  coordinates  [l8]   r,  9,   ^,   we  have 

(IX.l)    i^  .  ^  +  1,  11^  +  i  i  =  ^(p  +  (|  -c2-  |)T.+  o.u 

2irr   sm  0 

for  X-,  >  0   (for  x-,  <  0,  6(0)  is  replaced  by  6(0-Tr)  in  the 
source  term).   Here  l^  =  /l^  +  X^  +  i^  =   /|^  +  "X^  and 
Iq  =  tan"  "X/l-i^. 

Let  F(g,^  ,^0,r,0)  =  0  be  a  solution  of  (IX. 1).   Then 

II  =  _  ||/|l^  and  II  =  -  If/I^  and  equation  (IX. 1)  becomes 


The  transformation  of  a  multi-dimensional  6-function  from 
Cartesian  to  curvilinear  coordinates  is  discussed  in  [l8], 
pp.  291-293. 
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:ix.2)     fr^.^i. 


G 


{r,eA)  -l+o{i)  ^{^-^)^{^) 


L 


27rr  sin  0 


|£=  0 

Off 


J  og 


where  G(r,0,|)  ^  p  +  (|-  c^-  ^)T  +  c-u. 

Applying  the  method  of  characteristics  to  this  equation  for 
F,  we  have 


(IXO) 


dr  _  rde 


ag 


CO  ^ 


I  +  g(^)6(r-a,e) 
2irr   sin  6 


which  yields 


(IX. 4)  ||+-^=i 


CD 

L 


G 


,    a(g)5(r-a,0) 

2 
27rr   sin   9 


The  solution  Is 


r'/Ll 


(IX. 5)    g  = 


^  27rr  sin  6 


v/lI. 


dr 


J 


using  the  boundary  condition  at  oo,  namely,  g  =  0  for  r  — >  oo. 

Since  we  only  need  g  up  to  0(1/l),  we  expand  the  exponentials 
and  obtain 


(IX. 6) 


g  =  1^0)   ^   1  ~(1)  ^  _a^(,(^^0^^) 

^r 


a(e  ,^0) 
27r  sin  0 


2v     ^  .^r2 


^-i^        Llja   a^Ll^ 
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where,  in  G{r,9,i),   we  use  expressions  (VIII.2)  to  (VIII.6). 

The  function  g^  '{9,a,i)   describes  the  perturbed  velocity 
distribution  function  at  the  surface  of  the  sphere.   Granted 
that  the  object  is  sufficiently  small,  it  is  the  correct  first 
order  perturbation  and  takes  into  account  the  effect  of  all 
collisions  without  any  approximations. 

We  may  also  regard  g^  '^(9, a, I)  as  describing  the  distri- 
bution arising  from  a  modified  source  situated  at  the  center  of 
the  sphere  of  radius  a;  this  source  gives  rise  to  the  pertur- 
bation in  the  drag. 
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X.   Drag  on  a  Sphere  at  High  Mach  Number  and  for  Diffuse 
Reflection 


We  find  the  drag  by  calculating  the  momentum.  In  the 
direction  u  ,  Imparted  to  a  small  sphere  of  radius   a   sur- 
rounding the  source.   The  momentum  due  to  the  Impinging 
molecules  may  be  found  directly  from  g.   That  due  to  the 
reflected  stream  may  be  found  by  applying  the  appropriate 
boundary  conditions  on  g  at  the  surface  of  the  sphere. 

For  diffuse  reflection  we  have  already  calculated  N,  and 
N„  to  lowest  order  for  the  free  stream  case  (expressions 
(IV.4)).   We  now  find  N^  and  N^  to  0(1/l)  from 


r 


(X.l) 


N. 


(1) 


N 


(1) 


oo 

r 

00 

dl. 

de,e,  e  -=- 

0 

r 

di. 

00 

di. 

00 

r 

J 

0       —00     —00                     — 

00      - 

00 

-00 

0 

dl, 

^00 

dig 

.00      ^ 

00 

dl, 

dl^ 

J 

-i 

■00 

o  •^- 

00 

-00 

(1) 


(1) 


where   g^    ''    Is   the   term  of   0(i/l)    In  g. 


For  u     — y  00  we   have 


(X.2) 


where 


n(i)  -  0 


N 


(1) 


2   A(u    ,Q) 


2        2 
2TrL  r  cos   0 


A(u^,0) 


^^^^^    4°)r(|)  F(|.  1;  1;  -tan^e)   , 


(27r) 


3 


T 
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with 


n(°)  =  .^4^  Ve  u  cos  0  for  u  ^c» 


Thus,  the  reflected  stream  has  the  form  NX  '^e      where  y  =  T  /T. 
(see  Section  V) • 

The  drag  due  to  the  impinging  stream  in  the  u  direction  is 


r      r 


O  00  00 


(X.3)        D.  =- 


dS 


sj  \J 


dt 


di. 


di^I^(^^cos   9 


—  00  —00  —00 


+  ^2^1^  0)g(a,0,^^+ u^cos   6,   I^  +  u^sin  9,   \^) 


Here  dS  is  any  element  of  surface  on  the  sphere  whose  outward 
drawn  normal  is  in  the  positive  x.,  direction.   Because  the  flow 
is  symmetrical  about  u  ,  we  may  always  rotate  the  coordinate 
system  about  this  direction  so  that  u  lies  in  the  x-,Xp  plane  (or 
X, x^  plane).   We  note  that  in  calculating  the  drag  due  to  the 
impinging  molecules  we  need  only  consider  those  hitting  the  front 
of  the  sphere  when  u  — >  oo. 

The  drag  due  to  the  reflected  stream  is 


„00     ^00       ,00 


(X.4) 


D   =  - 

r 


dS 


df 


dt 


J 


de3(N(°) 


O     —00     -00 


+  ^)^^')e  ^     ^     ^  ^^(^^cos  0  +  ^2^1"  ^) 


again,  only  for  the  front  of  the  sphere.   The  total  drag  is  then 
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D  =  D.  +  D^ 
1    r 


If  we  let  L  -^  00,  we  recover  precisely  the  free-molecular 
flow  drag.   For  L  finite,  the  calculation  of  the  drag  reduces 
to  an  Integration  over  0   between  the  limits   0  and  Tr/2. 
The  integration  over  ?  yields 

r 


(X.5)  < 


D   = 
r 


D. 


T       yd 


T 


5/2 


dS(N^°^  +  N^^^cos  0 


2-3 

2-2   a  2  ■^^  ^o 
-  Tra  u  +  -r  Y 
o   L  ' 


r 


ir/2 


2iT  /tt 


sm 


F(|,l;l;  -  tan^e)d0 


o 


and  n|;  '    is  given  by  (X.2). 

Using  expression  (6)  of  Appendix  V,  we  may  evaluate  the 
integral  in  D.  which  becomes  the  dominant  term  for  u 
The  final  expression  for  the  drag  is  then 


(X.6) 


D  =  -  Tra 


■u. 


a 
L 


2-3 

y  ^o 

8  /2  J 
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XI.   Evaluation  of  g(x),  the  Perturbed  Distribution  F-unctlon, 

for  a  Plate  at  any  Mach  Number  and  for  Specular  Reflection 


Using  the  same  method  as  for  a  sphere,  we  solve  equation 
(III. 8)  for  g  using  polar  coordinates  x^,  9,    and  R.   In  the 
case  of  the  plate,  we  are  interested  In  the  various  moments  at 
X-,  =  0  for  the  purpose  of  calculating  the  drag.   This  leads  to 
considerable  simplifications  In  the  expressions  for  the  moments 
and  will  enable  us  to  evaluate  the  drag  for  all  Mach  numbers. 

Using  expression  (4)  of  Appendix  V,  we  find  that  the 
contributing  terms  for  the  plate,  assuming  specular  reflection, 
are  the  following  (up  to  0(1/L)): 


p(0,R) 


l67r/?  Au 
(27r)V 
327r^Au, 


o 


T(0,R)  = ^      (XI. 1) 

3(2Tr)^R^ 

u^^^O,R)  =  0 


Figurt  4 


for  X,  =0  (I.e.,  In  the  plane  of  the  plate).    In  polar 
coordinates,  the  source  term  a(^)5(x)  becomes 

^  ^   ^   5(x-,)5(R-a) 


27rR 


We  do  not  need  u 
case. 


WT 


and  u 


(3) 


to  evaluate  the  drag  In  this 
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The    solution   Is 


(XI. 2)  g{l^,y.  ,    R) 


-R/Ly 


O) 


G{0,R,l-^,iC) 


+ 


a(|^,X)5(x^)6(R-a)  I     ^^^-y_ 


dR 


27rR 


In  order  to  eliminate  the  x-,  dependence  In  the  source  term 
(this  X-,  dependence  has  already  been  eliminated  In  the  term  In 
G  by  using  the  values  of  p ,  T  and  u  In  the  plane  of  the  plate), 
we  use  the  conditions  on  the  characteristics,  namely. 


dx 


1   dR 


'1 


Integrating  the  source  term  over  x-,  yields 

5(e-L,:^)5(R-a) 


(XI. 3) 


and 


a(l^,:^)6(R-a; 
SttrL 


6(x-|_)dx^ 


27rR?. 


g^°^|^,"2,R) 


;xi.4)    i(!^,X,R 


03 


Au 


+ 


lX  (2Tr)^R 
S(|^,j^)5(R-a) 
2TrRL 


-16t>^  +  (1  H2  -  |)  2|L 


2  1 


The  first  order  correction  to  the  drag  Is 


It  Is  easily  seen  that  the  drag  due  to  g^    '  {i^,)C ,    R)  Is 
precisely  the  free  molecular  flow  drag. 
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(XI. 5)     D 


(1) 


4Au 


(a.)^(.:)^/^Lr°"-i-^^'-^'- 


-i~< 


u 

r  o 


~2 


+  [k^(l+u^)  +  2k^] 


-p  A 


e  2^  dx  -  Su^Ck-^  +  Sk^) 


where 


k,  =  -  l67r  /tt 


k^  = 


32  TT 

3 


At  high  Mach  number,  the  total  drag  approaches 


(XI.6)    D  =  -Au2  .  a  16]^  ^3  . 


o  ■  L  (2^)2    o 


Here,  as  for  the  sphere,  collisions  between  particles  bring 
about  a  reduction  in  the  drag. 
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Appendix  I:   Evaluation  of 


,2 
4-n    -t    ,, 
t   e        dt 

t-z 


n  =   1,2,3,^. 


The   case   of  n  =  0  has  already  been  explained   in  the   text, 
All   cases  with  n  >   1  may  be   reduced  to  the   n  =  0  case. 

1)        Ii(z)    =        ^Vi^  = 


J 


(t-z )  +  z       -t     , .  ,    . 

^    (t-z)      ^        d^  z   =  x+ly 


J 

-00 


yielding 


(1.1) 


lj^)(z)    =  e   ^  /tt  I  1  /tt  z  +   e^    -  2zl(z)  y>0 


l|^^(z)    =  e    ^  /^  1-1  Kir  z  +   e^    -  2zl(z)        y<0 


where 


I(z) 


e^   dp 


2)        I^(z) 


/-°°  2  r°° 

^2    -t     ,,  /,2         2^    ^    ^2       ,2 

tedt   ^  (t    -   ^,)    +   ^      e^   dt 

(t-z)  (t-z) 


(t-z) (t+z)   +   z^   ^-t^ 


T^^^ 


e        dt 


Thus 


(l.ii) 


/ 


I^^)(z)    =  e"^  Xt   U/?  z^  +   ze"^    -  2z2l(z)l     y  >   0 
42)(z)    .  e~^  K^  1-1/^  z^  +   ze"^    -  2z2l(z)!    y  <   0 
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3) 


13(2) 


00  00 


Tt^FT 


2         2 

't-z   +    z)t         -t    ,, 

^ TT H e        dt 

(t-z) 


yielding 


I-^^^(z)  =X7r  e"^  |)/^  iz^+^   (l  +  2z^)  -  2z^l(z)      j       y>0 


(l.iii)    { 


l^^\z)=yw  e-^^\-yT  iz^+\ 


v^ 


(l+2z^)- 2z^l(z) 


y  <    0 


^) 


14(2)    = 


00  ^  00 

t   e        dt 


Tt^FT 


(t-  z+  z)t^      -t^^, 
(t-z)^        ^        ^t 


yielding 


/^ 


l(^)(z)  =3yif  ze~^  Ix?  iz^+%-   (1+  2z^)  -  2z^l(z) 


(l.iv)      < 


K 


2(  2r  ) 


IjJ^Mz)=/^  ze"^    \-Y^  iz^+     2 


y  >  0 


y  <  0 


Integrals  with  values   of  n  >  4   are  not   needed   in   our  calcula- 
tions and  are,    therefore,   not   tabulated. 

The   remaining   coefficients   of  the   matrix    (ill. 11)    are 


a 


12 


2 


r 


&-  3)Qd! 


/-   .    2 


2iL  k 


"    1-2)/?  i  +  y?  \.a     +  ae"'  +   2(2-a'^)l(a) 


""21 


^ 
3 


J 


(c^-  3)fid|  =  I  a^2 


a  = 


/2 


u  k-,- 
o   1 


•L  k 


k 


46    - 


( 


^13 


J 


~        ~  k-,  2  (       _  2 

c^fid^    = -T^  e~^    <  I/tt  a  +  6°^    -  2al(a) 

IL  kl 


^31  ~  ^13 


'l4 


J 


^        ~  kp  _   2r  2 

c   fid^   = ^  e  °^    {  I/tt  a  +   6°^    -  2al(a) 

ILlkl 


a.i,  T    —  a 


ifl  -  ^14 


r 


a 


15 


c^fidC 


^ —  e  "^    ^  i/tt-  a  +  e"    -  2al(a) 


IL  k 


^51   ~  ^15 


a 


22         5 


-a 


(C^-      3)^nd|      =  ®     ^  1^     ^/-     -•  :.     ^/^     .._2       .      w-      .       ^ 


6iL|k| 


5  /tt  1-  4  /tt  la""  +  /-ir  la 


+  e"    [-  J  a  +  a^]   +  2l(a)[-5  +  ^a^-  a^] 


a 


23  ~        3 


k.  e 


-a 


( 
c^('c^-   3)nd|   =    -  \  ^7^7-7^1  -  2/?  la  +  X?  la^ 


lL|k| 

2 

+  e°^    (-  I  +  a^)   +  2al(a)[2-a^] 


^32 


2 


~        ~2  ~  T 


'•24 


3 


,~2 


k^e 


-a 


Q-Ac    -  3)fide   -   -  4 


1   '^2 


3    .^  I,   ,2 
iL  k 


-  2  /tt  la  +  y?  la^ 


2 
+  e"    (-  I  +  a^)   +  2al(a)[2-a^] 
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r 


a 


^2 


a 


25  3 


-a 


2  /tt  la   +  K"""   let 


IL  k 


+   6°^    (-  I  +  a^)    +   2al(a)[2-a^] 


r 


a 


52   ~    "   3 


c^(c^-   3)fid|   =  I  a2^ 


a 


35 


~2         ~            p"°^^      I  ^^2    +    ^3)    ^^   ^         w-  2    2 

c,  fi  d^   =  ^ ^  {  —^ ^ +  Ktt  1   k^a^ 


'1 


IL  k 


3 


+  k^ae°^    -  21(a) 


r     2  2 

(ko  +  k^)  22 


^2         ^^3' 


+  k-,a 


a 


3^ 


k,k^  2 

~    ~    r^^£  12         -a 

^   2  ILlkp 


&1    .    ./^..2 


+  /tt  la     +  ae 


a 


+   I(a)[l-2a    ] 


^35 


~   ~   n^?  ^1^3        -a^  )      /il    .    ./=  --.2 

^   5  ILlkP  '  2 


+  /it  la     +  ae^ 


+   I(a)[l-2a^] 


a 


53 


3  1 


=  a 


35 


r 


^^3 


CgC-j^fid? 


=  a 


35 


r 


'43 


c2c-|^ndC 


=  a 


3^ 


-  i^8   - 


'-44 


+ 


\^l 


(k^   +   k2)     '     (k^   +   k2)|k|2j 


2  2  2 

kp        _       p      k„  a  e°^ 

+  — %  Ktt  la    +  — — 5 21(a) 


k^k| 


^3' 


\C 


+ 


+ 


2,  2  -, 
a  kg 


2 (kg   +   k|)|k|^    '    2(k^   +   k|) 


r 


'■45 


2    p 


-a 


Wit  1 


k^kgk^ 


kgk^ 


iL|k|     [    2       L(k2+k|)|k|^      (kgH-  k|)  _ 


2kpk^  p        kpk^ae 

+  -^-4  /tt  la^  +  -^^-p —  +   I  (a ) 


k 


2k„k^a 
2    3 

ik|2 


LkL 


kgk^ 


k-^kgk-j 


.(k^  +  k|)         (k^  +  k|)lkl2 


^5^ 


c^CpfidC   =  a^i^p- 


a 


55 


J 


c| fid^   =  — 


-a 


)y^i 


ILjkl 
2 


k' 


2- 
^3' 


+ 


k^k^ 


L(k^  +  k^)     (k^  +  k2)|k|2  J 


y-w  ik^     2  S  a^ 

+  ^  a'^  +  — %  ae"    -  21  a 


1  2    2 
k^a 

Ik|2 


1  2,  2 

1  :> 


^y 


+ 


_2(k^   +   k|)lk|^     '    2(k^   +    k|) 


The  values   of  these   coefficients  at  high  Mach  numbers  are   listed 
In  Appendix  II. 


In  calculating  the   a.  .  we  use   the   relations 
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I     k         kp       k^ 
^    Ikl    1    Ikl   2    k   3 


%l     +     k|  /k^  +  k| 


I   I   I 


expressing  the  new  coordinates  i-^,    i^,    i^     in  terms  of  the  old, 
The  Inverse  relations  are 


2  .  ,  2 


k       fe^^;  .' 

^^"Ud"^i"   Ikl   '^ 


kp   ,       1^^     -       ^1^2     ^  ' 
1^1      j/k^  +  k|       pr+4  li^l 

'   1^1  /i^^T^      /i^fT^lkl 
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Appendix  II 


(a)   Val-ues  of  the  coefficients  a.  .  and  the  source  terms  a(k) 


when  u  — >  oo  . 
^  o 


The   leading  terms  to  0(l/L   )    are 
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^  L|kp  [  2  ^^   L|k| 


/2    T.lk 


Llkl         /2 


lL|kr  V  y2    |k| 


-(2).,        ^^1^24°^ 
-^  L  kP 


2         L|k|    1        lL|k|      \  /2    Ikl 


^f)(k; 


o^{\^) 


Ak^k^N^ 


o) 


L  k 


3 


.  ^  .  /I^l  .  ^^^^ 


+ 


1   - 


y2Y  N^°^Ak-j_ 
31L^|kp 


2         L  k  IL  k 


~2 
-,       /^s      Ak-,  Au^ 

2    ^ T         ,.  i,_i2    ^ 


)/2    Ik! 


3 


1     ^ 


lL|k|^        lL|k|    V3/2        6    |kp 


Au 


o  1 


k: 


3L^|kpV  X2 


k 


for  diffuse  reflection. 


(b)   Source  Terms  for  any  Mach  N-umber  and  Specular  Reflection 


1        L|ki       \    lk| 


+  u 


k. 


^(l)rv^   2Ae""  [  g^/  ^1    ^~   '^1  ^   o 


IL 

1  2  , 

k  k 

+  — ^  a    -   2a  — ^ 


k 


k 


/         1    o  Is 

(a  = — ) 

/2  |k| 


(k^  +  k|) 


2  k 


3 


k 


3 


k 


1(a) 


2Ak. 


1(2)  (ic)  =^^2^e- 


2  (  k- 


3(3'(k) 


lL|k 
2Ak 


3   3-a 


k 
2  r  k 


[ae"  +  (l-2a  )l(a)]  +  I/tt  u^a 


IL  k 


3 


^1,^2  _    ^ 

— ^  [ae"  +  (l-2a  )l(a)]  +  I/tt  u  a? 
Ikl  °  ) 
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Appendix  III:  Some  Properties  of  the  Modified  Bessel  Functions 


The  modified  Bessel  functions  of  the  first  and  second 
kinds,   1^(2)  and  K  (z)j are  defined  as  solutions  of  the 
equation 


dz 


We  have  I,,(z)  =  e 


•V7ri/2 


J  (iz)   which  is  a  real  function 


V 


of   z,    and     K^(z)    =  |  e^^+^^^^Z^   h(1)(iz) 
For      V     not   an   Integer,    we   have 


I      (z)    -   I    (z) 

K  (z)  =  2  ^ : 

Sm     VTT 


For      0  <  X  <   1 


K    (z)  «  ^n  4- 
Y 

Kjz)^  \   r(n)(|)"        n  =   1,2, 


Asymptotically, 


(7  =  Euler's   constant) 


K^(z)--.v/^  e 
^  /2z 


-z 


For  further  properties   of   these   fvinctlons   see    [I9]    and   [20] 
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K  (z)  exists  for  all  finite  values  of  v  and  Is  monotonlcally 
Increasing  with  v  for  a  given  z.   For  all  values  of  v,  however. 


Kv(0) 


The  K  (z)  satisfy  the  following  relations; 


1)  s_i(^)  -  Vi(-)  =  -l^s(-) 


2)        \.i(^)    +   VK^(2)  =  -2K;(z) 


3)        zK^(z)  +  vK^(z)  =  -zK^_^(z) 

m  (  K  (z)  ")         K  ,  (z) 


v+m 
z    j         z 


3)    (^ri^h<-) 

6)  k;(z)  =  -Kj{z) 


7)         K_^(z)  =  K^(z) 

dx-. 
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Appendix  IV:   List  of  Fourier  Transforms 

All  the  following  transforms  are  listed  In  Magnus,  W. , 
Oberhettlnger,  F.,  and  Tricoml,  F.  G.,  Tables  of  Integral 
Transforms,  California  Institute  of  Technology,  Bateman 
Manuscript  Project,  A.  Erdelyl,  editor  (McGraw-Hill  Book  Co, 
195^). 


1) 


2) 


3) 


M 


5) 


2^- 


x2(a''  +  x^)  2  (xy)2  J^(xy)dx  =  y  ^e 


-2^ -ay 


Re  a  >  0 
y  >  0 


Vol.11,  p. 7,  No. 4 


-ax 


X 


sln(xy)dx  -   tan   (a  y) 


o 


Vol.1,  p. 8,  No. 11 


^°g(^y)   dx  =  I  a-^  e-^y 


Vol.1,  p. 8,  No.18 


^  /   \  -ax,     /  2   2  N  --2 
J^(xy)e   dx  =  (y  +  a  ) 


o 


Re  a  >  0 
y  >  0 

Re  a  >  0 
y  >  0 

Re  a  >  0 
y  >  0 


Vol. II,  p. 9,  No.18 


sln(xy)  — 2 —  2  V2  ^^  "  Y  KQ(ay) 
(a  +  X  )^' 


Re  a  >  0 
y  >  0 


Vol.1,  p. 66,  No. 27 
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6) 


4 


.2  2,-i 


xJ  (xy)KOx)dx  =(p^  +  y  +  2p  y  ) 


Re  P  >  0 
y  >  0 


7) 


8) 


9) 


10) 


-V+i,  2 


V-ii,  M,+^ 
(x^  +  a^)^-"j^(xy)(xy)*dx  =  l^\-^\^+i)    ^v-^l^^^^ 


2,^1-1. 


o 


Re  a  >  0     y  >  0 

-1  <  Re  V  <  2Re  \i  +   3/2 

Vol.11,  p. 24,  No. 20 

^^   K,(ax)cos(xy)dx  =  >^  (2a)-^'r(tli  +  i)(y^  +  a2)+^-* 


Re  p.  >  -|  if  upper  signs  are  used 

Re  i-L  <  -g-  If  lower  signs  are  used 

Re  a  >  0      y  >  0 

Vol.1,  p. 49,  No. 41 


''2'H^nl^  ^^  =    (*  |)"[r(n+i)]-^   Kjay) 


(a"+  x^)^ 


Rea>0,      Ren>-i,      y>0 
Vol.1,   p. 11,    No. 7 


x       K    (ax)sin(xy)dx 
I-'- 


-  yr(i^i-i>^  +  i)r(i-*^  -M  F  (^±^,  ^4^;  |;-^) 

Re(At^i)   <2,      Rea>0,      y>0 
Vol.1,   p. 106,    No. 50 
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11) 


r  cos(xy)  =  (-^^L^    ^^    ^  [yX(^y)l 


o 


(x    +  a    )      ^ 


2V(v+i)       d/ 


Rea>0,      0<m<Rev+^,      u>0 
Vol.1,    P. 14,    No. 28 


12) 


0 

J 


x'^-2(xy)2j^(xy)K^(ax)dx 


13) 


_  r[i(v-A+^+i)]r[^(v->--^+i)]  p  /v-x+Eti  v-^-^^+i.  v+H .  y?.^ 


2'+v-^+ir(v+i)y-'-^ 


r  2m+l 

X 


Re   a  >   0,      y  >   0,   Re(v-A+l)>  |Re(iJ.) 
Vol.11,   p. 63,   No. 4 


/    -,  sin+l./—        -,2m+l 
(-1)        yir       6_ 


o 


(a    +x    )      2  2   a    I     (n+2)      dy 


-2  <  2m  <   2n,      Re   a  >   0,      y  >   0 
Vol.1,   p. 67,    No. 37 
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Appendix  V:   Some  Properties  of  the  Hyper geometric  Functions 
F(a,b;c;z) 

The  hypergeometrlc  functions  are  solutions  of  the  differential 
equation 

2 

(1)  z(l-z)  ^+  [c  -  (a+b+l)z]  ll  -  abu  =  0  . 

dz 

One  of  the  solutions  of  (l)  may  be  represented  by  a  series  which 
Is  regular  at  z  =  0.   This  "hypergeometrlc  series"  Is 

(2)  u   Ffa  b-c-z)    n  I  ab  z    a(a+l)b(b+l)  z 

When   c  =  -n  (n  =  0 ,1 ,2  ,3 ,  .  .  .)  ,    the  denominator  of  the  (n  +  2) 
term  and  the  following  all  vanish  and  the  series  Is  undefined. 
However,  the  following  representation  Is  valid: 

c  =  -n    I  '(c) 

a(a+l) . . . (a+n)b(b+l) . . . (b+n)F(a+n+l,  b+n+1;  n+2;  z)  ^n+1 

(n+1) I 

The  hypergeometrlc  series  converges  on  the  unit  circle  \z\    =1 
In  the  following  cases: 

a)     1  >  Re(a+b-c)  >  0     convergence  on  the  entire  unit  circle 

except  for  z  =  1 
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b)    Re(a+b-c)  <  0     absolute  convergence  on  the  entire  unit 

circle  Including  z  =  1. 

If  Re(a+b-c)  >  1  the  series  diverges  on  the  entire  unit  circle. 
In  all  these  cases  we  must  assume  that  a,  b,  c  are  different 
from  zero  and,  furthermore,  that  they  are  not  negative  Integers. 

A  transformation  formula  for  the  hypergeometrlc  function 
which  is  frequently  used  In  this  report  Is 

(4)    F(a,b;c;z)  -  r(c)r(b-a)  (_z)-a  p^^,  1-c+a;  1-b+a;  ^) 

r(b)r(c-a) 
r(a)r(c-b) 

We  also  use 


(5)    11m   P(a,b;c;z)= 
c  -^  -n 

=  a(a+l)...(a+n)b(b+l)...(b+n)z"+^  p(a+n+l,b+n+l;n+2 ;z ) 

(n+1)! 

-a 


P(a,b;c;z)  =  (l-z)    F(a,  c-b;  c;  ^^j) 

(6)  [ 

=  (l-z)  ^  F(b,  c-a;  c;  -^) 
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